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Abstract 

^^ In this paper we present a modification of a technique by Chiba and Nishizeki [Chiba 

Z^ . and Nishizeki: Arboricity and Subgraph Listing Algorithms, SIAM J. Comput. 14(1), 

pp. 210-223 (1985)]. Based on it, we design a data structure suitable for dynamic graph 
algorithms. We employ the data structure to formulate new algorithms for several prob- 
lems, including counting subgraphs of four vertices, recognition of diamond-free graphs, 
cop-win graphs and strongly chordal graphs, among others. We improve the time com- 
pel ' plexity for graphs with low arboricity or /i-index. 

Keywords: data structures, dynamic algorithms, arboricity, ft.-index, cop- win graphs, 
-. diamond-free graphs, strongly chordal graphs. 
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c/3 ! 1 Introduction 

We describe a variation of a technique by Chiba and Nishizeki [3] , leading to a data structure 
for graph algorithmic problems, called the h-graph data structure. It supports operations of 
insertion and removal of vertices, as well as insertion and removal of edges. Although the 
data structure can be used for general purpose, it is particularly suitable for applications in 

^^ I dynamic graph algorithms. 

As an application of this data structure, we describe new algorithms for several graph 
problems, as listing cliques; counting subgraphs of size 4; recognition of diamond-free graphs; 

f^ I finding simple, simplicial and dominated vertices; recognition of cop-win graphs; recognition 

of strongly chordal graphs. We remark that no previous such dynamic algorithms exist so 
far in the literature. On the other hand, in some cases, there is also an improvement in time 

k> I complexity (relative to static graph algorithms) for graphs of low arboricity or /i-index. 

^ ' A dynamic data structure designed for graphs with low /i-index has been first defined by 

Eppstein and Spiro [X\- This data structure keeps, for each graph G with /i-index h, the set of 
vertices with degree at least h, and a dictionary that indicates the number of two-edges paths 
between any pair of vertices, for all those vertices at distance 2. The total size of the data 
structure is 0{mh) bits. Using this information, the authors show how to maintain the family 
of triangles of G in 0(h) randomized amortized time while edges are inserted or removed. 
The authors also show how to keep other statistics of G with this data structure. Our /i-graph 
data structure follows a different approach. First, we store no more than the adjacency lists of 
G in a special format, using 0{n + m) bits. Second, we do not compute the /i-index of G. The 
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advantage of our data structure is that we can use to maintain the family of triangles of G 
in 0{dh) deterministic worst case time per vertex insertion or removal, where d is the degree 
of the vertex. Furthermore, the time required by this algorithm when it is applied to all the 
vertices of the graph, so as to compute the family of triangles of G, is 0{am), where a < h 
is the arboricity of G. The disadvantage is that we cannot longer maintain the triangles 
as efficiently as Eppstein and Spiro when edge operations are allowed. So, though both 
data structures have some similarities in their inceptions, they are better suited for different 
applications. In particular, our /i-graph data structure allows the "efficient" examination of 
the subgraph of G induced by the neighborhood of an inserted or removed vertex. 

One of the similarities between the /i-graph data structure and the data structure by 
Eppstein and Spiro, is that both differentiate between low and high degree vertices. The 
technique of handling differently vertices of high and low degree has been first employed by 
Alon et al. [1], and since then many other works made use of this classification (e.g. [71llU|). 
However, we use a local classification on each vertex. So, some vertices can be considered as 
both high and low depending on the local classification of each of its neighbors. 

The paper is organized as follows. In the next section we introduce the notation and 
terminology employed. In Section [3] we discuss the technique by Chiba and Nishizeki, and its 
variation for dynamic graphs. The /i-graph data structure is described in Section HI together 
with the operations that it supports. Finally, in Section [5l we show the applicability of the 
h-graph data structure by solving the problems listed above. 

2 Preliminaries 

In this paper we work with undirected simple graphs. Let G be a graph with vertex set V{G) 
and edge set E{G), and call n = \V{G)\ and m = \E{G)\. Write vw to denote the edge of 
G formed by vertices VjW £ V{G). For v £ V{G), represent by Ng{v) the subset of vertices 
adjacent to v, and let A''G'[f] = Ng{v) U {v}. The set Ng{v) is called the neighborhood of v, 
while A'^q[t;] is the closed neighborhood of v. The edge-neighborhood ofv, denoted by Nq(v), is 
the set of edges whose both endpoints are adjacent to v. Similarly, the neighborhood Ng{vw) 
of an edge vw is the set of vertices that are adjacent to both v and w. All the vertices in 
Ng{vw) are said to be edge-adjacent to vw. The degree of v is dciy) = \N{v)\, the degree of 
vw is dcivw) = \Ng{vw)\, and the edge-degree of v is d'Q{v) = \Nq{v)\. When there is no 
ambiguity, we may omit the subscripts from N and d. 

For W C V{G), denote by G\W] the subgraph of G induced by W, and write Eg{W) to 
represent £^(G[VF]). As before, we omit the subscript when there is no ambiguity about G. 
A clique is a set of pairwise adjacent vertices. We also use the term clique to refer to the 
corresponding induced subgraph. The clique of size k is represented by K^, and the graph 
K^ is called a triangle. We shall denote by 0{n^) the time required for the multiplication of 
two n X n matrices. Up to this date, the best bounds on n^ are v? < n^ < n?-^"^^ [5]. The 
arboricity a{G) of G is the minimum number of edge-disjoint spanning forests into which 
G can be decomposed. The h-index h{G) of G is the maximum h such that G contains h 
vertices of degree at least h. It is not hard to see that 

- < — — < a{G) < h{G) < V2^ 
2 n — 1 

for every graph G, where 6 is the minimum among the degrees of the vertices of G. 



For each vertex v of a graph G, define N{v,i) = {w £ N{v) \ d{w) = i}, i.e., N{v,i) 
is the set of neighbors of v with degree i. Denote by L{v) the set of neighbors of v of 
degree at most d{v) — 1, and H{v) the set of neighbors of v of degree at least d{v), i.e., 
L{v) = N{v, 1) U . . . U N{v, d{v) - 1), and H{v) = N{v, d{v)), ..., N{v, n - 1). We use i{v) 
and h{v) to respectively denote |^(w)| and \H{v)\. Observe that v can have at most h{G) 
vertices of degree at least d{v) + 1, thus h{v) < h(G), and the number of nonempty sets in 
the family N{v, 1), . . . , N{v, d{v)) is at most 2h{G). 

3 Revisiting the Approach by Chiba and Nishizeki 

In [3], Chiba and Nishizeki devised a new method for listing all the triangles of a graph, based 
on the following lemma. 

Lemma 1 ([3]). For every graph G, 

y Tii\n{d{v) , d{w)^ < 2a{G)m. 

vweE{G) 

The idea of the algorithm is simple; for each v G V{G), find those z G N{w) that are 
adjacent to v, for every w £ L{v)DN{v, d{v)). (The original algorithm by Chiba and Nishizeki 
is slightly different, so as to list each triangle once. In particular, it requires the vertices of 
V{G) to be ordered by degree.) By Lemma [U such algorithm takes 0{n + a{G)m) time. 

At each iteration, the algorithm by Chiba and Nishizeki finds only those triangles v, w, z 
such that d{v) > nim{d{w) , d{z)} . Suppose, the aim is to find all those triangles containing 
some vertex v. For instance, suppose we need to dynamically maintain all the triangles of 
G, while vertices are inserted into G, and v has been recently inserted. In this case we can 
find those z £ H{w) that are adjacent to v, for every w G N{v). The total time required 
for computing all the triangles with this algorithm while vertices are dynamically inserted is 
again 0{n + a{G)m), according to the next lemmas (see Section [3] for the implementation 
details). 

Lemma 2. Let ei,. . . ,em be an ordering of E{G) for a graph G, and call Cj = ViWi. Denote 
by hi{v) the value of h{v) in the subgraph of G that contains the edges ei, . . . ,ej, for every 
1 < i < m. Then, 



J2hi{vi)<MG)m. 



i=l 

Proof. Call di (v) to the degree of v in the subgraph of G that contains only the edges ei , . . . , Cj . 
We define two values, fi{vw) and F{vw), for every pair of vertices v and w that are useful 
for decomposing the values of /ij. Specifically, for 1 < i < m, let 

Jl iiv = Vi, w £ NG{v)\{wi} and di{w) >di{vi) 
fi[vw) = < 

I otherwise 

and F{vw) = J27^i fiiy^)- With this definition, 



i=l i=l vweE{G) vweE{G) 



We now prove that F(wtt;) < 2inm{dG{v),dGiw)}. For this, suppose that dG(w^) — '^G(^) = 
k, for some A; > 0. By definition, if Cj 7^ vw is inserted before vw, then fi{vw) = 0. 
Similarly, if ej 7^ vw is one of the last k — 1 edges inserted among those incident to w, 
then di{w) > dciv) > di{v), thus fi{vw) = 0. Therefore, at most dciv) edges ej incident to 
w are such that fi{vw) = 1, which implies that F{vw) < Idciv) as desired. Consequently, by 
Lemma [H 

m 

Y.^ii^'i)^ Yl ^i^^)< Yj 2min{dGW,dGM}<4a(G)m. 

i=l vweE{G) vwGE{G) 

D 
Lemma 3. For every graph G, > h{v) < Aa(G)m. 

vweE{G) 

Proof. The proof is similar to the one of Lemma [2j Just replace di with d in the definition of 
fi and follow the same proof. D 

One of the properties about this new algorithm for computing the triangles, is that we 
can apply it to each vertex v so as to find all the triangles that contain v in 0{d{v)h{G)) 
time. In Section 15.21 we also show how to count the number of some induced graphs on four 
vertices that contain v. 

The next corollaries are also relevant for the /i-graph data structure to be presented. 

Corollary 4. Let vi,...,Vn be an ordering ofV{G), for a graph G. Denote by hi[v) the 
value of h{v) in the subgraph of G induced by vi, . . . , Vi, for every 1 < i < n. Then, 

n 

Y Y ^^(^) - MG)m. 

j=i weN(vi) 
Corollary 5. For every graph G, > > h{vu) < 8a{G)m. 

v&V{G)w&N{v) 

4 The /i- Graph Data Structure 

Li this section we present a new data structure, called the h-graph data structure, that is well 
suited for graphs with low /i-index or low arboricity. The goal is to make it easier to traverse 
the set H{v) in 0{h{v)) time, for every vertex v, while vertices and edges are dynamically 
inserted to and removed from the graph. Then, this data structure can be used to solve many 
problems by using the technique implied by Lemmas [2] and [3] and Corollaries H] and \S\ as the 
maintenance of triangles described in the previous section. 

For each vertex v, the ^-graph data structure stores an object with the following data: 

• The degree d{v) of v, 

• A doubly linked list J\f{v) containing one object N{v, i) for each of the nonempty sets of 
the family {N{v, 1), . . . ,N{v,d{v) - 1)} (1 < i < d{v)). The members of the list J\f{v) 
are ordered so that N(f,i) appears before N{v,j) for 1 < i < j < d{v) — 1. 



• A set H(w) representing H{v). 

• A pointer data(t;), referencing an object that contains the data of v. 

For simphcity, we refer to H(f) as N(f , d{v) + z), and to H{v) as N{v, d{v) + i), for every i > 0. 
That is, for every 1 < i < n — 1, the set N{v,i) contains the neighbors of v with degree i, 
while N(f , i) is the object of the data structure that contains the neighbors of v with degree 
i. In the data structure, each }i{v,i) (1 < i < n — 1) is stored as a doubly linked list that 
contains one object for each w G N{v,i), with the following data associated to the object 
representing w: 

• a pointer Tpos{v,w), referencing the object of }i{w,d{v)) that represents v, 

• a pointer list(f,tf), referencing the list N{w,d{v)), 

• a pointer node{v,w), referencing the object that represents w, and 

• a pointer da.ta.{v,w), referencing an object that contains the data of vw. Note: both 
da.ta.{v,w) and da.ta.{w,v) reference the same object. 

Observe that the space required for the /i-graph is 0{n + m) bits. The /i-graph data 
structure provides the set of basic operations described below. Inside the parenthesis we 
show one value (x) or two values (x, y). The value x is the time required when the operation 
is applied once, while y is the time required when the operation is applied once for each vertex 
or edge. (We assume m > n.) For instance, the insertion of one vertex takes 0{dh) time, 
while the iterative insertion of the n vertices takes 0{am). Here d = d{v), h = h{G) and 
Q = a{G). 

• vertex_insert: inserts a new vertex v into G with a specified neighborhood N(v) 
{0{dh), 0{am)). 

• vertexjremove: removes a vertex v from G {0{dh), 0{am)). 

• edge_insert: inserts a new edge vw into G {0{h), 0{am)). 

• edge_remove: removes an edge vw from G {0{h), 0{am)). 

• adjacent: queries if two vertices v and w are adjacent {0{h)). 

• H: returns the set H{v) for a vertex v (0(1)). 

• N': returns the set N'{v) for a vertex v {0{dh), 0{am)) 

• G[N()]: returns the adjacency lists of the graph G[N{v)] for a vertex v {0{dh), 0{am)). 

The implementation of H is trivial. We discuss the other operations below. In the following 
sections we show several algorithms that work on the /i-graph data structure. 



The insertion of vertices and edges. The algorithm for inserting a new edge vw into G 
is straightforward. In a first phase, update the families J\f{z) for every z G N[v] U A''[t(;]. For 
this, create the set }^{v,dGiv)), move the vertices with degree dciv) from E{v) to ^{v,dG{v)), 
and move v from N(z, dciv)) to N(2;, dciv) + 1), for every z £ H{v). Next, apply the analogous 
operations for w. The second phase is to actually insert the edge vw. For this, insert v at 
the end of N(i(;, dciv) + 1) and w at the end of N(f , dciw) + 1), update the values of d{v) and 
d(u)), and create the pointers pos, list, node, and data for vw. 

Discuss the time complexity of the above algorithm. For the first phase, apply Algorithm[T] 
twice, once for v and once for w. Recall that this algorithm is applied before incrementing d 
for V and w, thus d(z) is the degree of z before the insertion of vw. Note that each iteration 
of the main loop can be implemented so as to run in 0(1) time, by using the pointers in the 
data structure. Thus, the update of AA(z), for every z G N[v\ U A^[tf], takes 0{h{v) + h{w)) 
time. 

Algorithm 1. Update oi N{z) for every z G N[v]. 

1. Insert a new empty set N(v,d(f)) at the end oi J\f{v). 

2. For each z G H(w): 

3. If d(2;) = d(w), then move z from E{v) to Yi{v,d{v)). 

4. Move V from N(z, d{v)) to N(z, d{v) + 1). If N(z, d{v)) = 0, then delete N(z, d{v)). 

5. liY\{v,d{v)) = 0, then delete Yi{v,d{v)). 

For the second phase, traverse the family J\f{w) until the first set N(w, d) with d > dciy) 
is reached. Then, create the set N = N(tt;, dciv) + 1), if d > d{v) + 1, and insert v into N. Next, 
traverse M^v) so as to find the set that must contain w, and insert w. Recall that there are 
at most 2h{G) sets inside each of J\f{v) and J\f{w). So, the time required by these steps is 
0(min{d(i;), d(u;), /i(G)}), while the creation of the pointers and the increase of d take 0(1) 
time. Therefore, the insertion of vw requires time 

O {mm{d{v),d{w), h{G)} + h{v) + h{w)) = 0{h{G)) 

The insertion of a new vertex v is simple. First, insert v as an isolated vertex, and then 
add the edges vw, for each w G N{y). The time required is 0(1 + d{v)h{G)). If we use 
the above algorithm for building G from scratch, then the total time is 0{n + a{G)m), by 
Lemma [2j 

Corollary 6. The time required for inserting the vertices and edges of a graph G, one at a 
time in no particular order, into an initially empty h-graph data structure is 0{n + a{G)m) . 

Proof. The insertion of the n vertices takes 0{n) time. Let ei, . . . , e^ be an ordering of E{G), 
and call di{v) and hi{v) to the values of d{v) and h{v) in the graph prior the insertion of Cj. 
By the analysis above, the time required for the insertion of Cj = vw is 

0(min{dj(u), di{w)} + hi{v) + hi{w)). 

Thus, by Lemmas [1] and [21 the total time required for the insertion of all the edges is 

Y, 0{mm{di{v),di{w)])+ Y. 0{hi{v) + h,{w)) = 0{a{G)m) 

vweE(G) vweE{G) 

D 



The removal of vertices and edges. For the removal of an edge vw, we should undo 
the insertion process. But this time, we first undo the second phase, and then we undo the 
first phase. To undo the second phase we need to physically remove the edge vw. Suppose 
that dciv) < dciw), i.e., w G H{v). Traverse E{v) so as to locate and remove the object 
that represents w in H(w). By using the pos and list pointers, we can easily remove v from 
N{w,dc{v)) in 0(1) time. Since v has h{v) neighbors in H{v), this phase takes 0{h{v)) = 
0{h{G)) time. To undo the first phase, we need to update the families M{z) for every 
z G Ng{v) U Ng{w). For this, we move v from N{z,dG{v)) to N(z,dG{v) — 1) for every 
z £ Hg{v), and then we remove ]^{v,dG{v) — 1) from J\f{v) so as to append it to H(r;). Note 
that in this step we should update all the list pointers of the vertices in N(f , dG{v) — 1) so 
as to reference H(u). Next, we should apply the same operations for w. This phase is rather 
similar to the one for edge insertion and it also takes time 

0{hG\{vw}{v) + hG\{vn,}{w)) = 0{h{G)). 

As before, we can remove a vertex v by removing all its incident edges first, and the object 
representing v later. The time required by this operation is again 0(1 + d{v)h{G)). Finally, 
if we use the above algorithm for decomposing G, then the time required for removing all the 
edges is 0{n + a(G)m), by Lemma[2j 

Corollary 7. The time required for removing the vertices and edges of a graph G, one at a 
time and in no particular order, from an h-graph data structure is 0{n + a{G)m). 

Adjacency query. Querying whether two vertices v and vu are adjacent takes 0{h{v) + 
h[w)) = 0{h{G)) time. Simply traverse the set H{z) for z G {v,w} with minimum degree. 

Edge-neighborhoods and their induced subgraphs. To compute the set N'{v) for a 

vertex f, first mark each z G N{v) with 1. Following, traverse each z G H{w) for every 
w G N[v) and, if z is marked with 1, then insert it into N'{v) and mark it with 2. Since 
each w is traversed in 0{h{w)) = 0{h{G)) time, this algorithm takes 0{d{v)h[G)) time. 
Furthermore, the time required to find {N'{v)}y^yiG)i by applying this algorithm to all the 
vertices in the graph, is 0{a{G)m), by Corollary [5l 

Clearly, the subgraph of G induced by N{v) is just the graph whose vertex set is N{y) 
and whose edge set is N'{y). Thus, the graph G[A^(v)], implemented with adjacency lists, 
can be computed in 0{d{v)h{G)) time, while the family {G[N {v)]} ^^y (^g) can be computed 
in 0{a{G)m) time. 

5 Applications 

In this section we list several problems that can be improved by using the technique implicit 
in Lemma [2] and Corollary HI These algorithms build upon the /i-graph data structure, and 
serve as examples of its applicability. One of the most appealing aspects of some of these 
algorithms is that they are simple to obtain from the problems' definitions. 

5.1 Listing the Cliques of a Vertex 

In [3] , Chiba and Nishizeki devised an algorithm for listing all the K}^ 's of a graph, for a given 
fc G N. Conceptually, their algorithm computes all the K^^is in the subgraph of G induced 



by N{v), and it outputs v plus these cliques, for every v £ G. (Again, some ordering of the 
vertices is required so as to avoid repetitions.) 

This algorithm can be translated to find all the K^s that contain a given vertex v. First, 
compute G' = G[N(v)] in 0{d{v)h[G)) time, as in Section HI Then, apply the algorithm by 
Chiba and Nishizeki that lists the Kk-is of G' in 0{\V{G')\ + {k- l)a{G'f-^\E{G')\) time. 
Note that 

\E{G')\< Y^ mm{d{v),d{w)} < 

weN(v) 

< Yl ^H+ Yl d{v) = 0{d{v)h{G)), 

'wGN{v) w£N(v) 

d{w)<h{G) d{w)>h(G) 

thus the time required to find all the -ftTfc's that contain v is 

O (d{v)h{G) + \V{G')\ + ka{G')^^^\E{G')\^ = O (^kd{v)h{G)a{G)''-^^ . 

Also, the total time required for listing all the cliques, by iteratively executing the above 
algorithm for each vertex v (and then removing v), is 



O [n + a{G)m + ka{G)'''^ Y Y ^^Hd{v) , d{w)} \ = O (n + ka{G)'' m 

\ v&V{G)weN{v) / 

matching the time complexity of the algorithm by Chiba and Nishizeki. 

5.2 The 4-Subgraph Counting Problem 

The 4-subgraph counting problem is the problem of counting how many copies of some graph 
H on four vertices appear as induced subgraphs of a given graph G. The connected graphs 
on four vertices are six: the complete graph K4, the diamond iC4 \ {e} for e € E{K4), the 
square C4, the path P4, the paw P3U Ki, and the claw K3 U Ki. The disconnected graphs on 
four vertices are five: K4, K4 \ {e}, C4, P3 U Xi, and K^ U Ki. In [TU], Kloks et al. showed a 
system of linear equations that solves the 4-subgraph counting problem for connected graphs. 
Specifically, Kloks et al. proved the following theorem. 

Theorem 8 ([10]). Let H he a connected graph on four vertices such that there is an 0{t{G)) 
time algorithm counting the number of induced H 's in a graph G. Then, there is an 0{n^ + 
t{G)) time algorithm counting the number of induced H's of G for all connected graphs H on 
four vertices. 

Since the number of K4S can be computed in either 0{n + m^'^^'^'''^) = 0{n + m^'^^) or 
0{n + a{G)'^m) time [3l [10], solving the 4-subgraph counting problem for connected graphs 
takes 0{n'^ + minjm,^'^^, a{G)'^m}) time. This result can be improved with a new system of 
linear equations that solves the 4-subgraph counting problem in 0{n + a{G)m + t{G)) time, 
even for disconnected graphs. By using the algorithm in [3j for counting the number of -ftr4's, 
an 0{n + min{m^'^^, a(G)^?n,}) time algorithm is obtained. The system of linear equations 
appears in the proof of the next theorem. 



Theorem 9. Let H he a graph on four vertices such that there is an 0{t{G)) time algorithm 
counting the number of induced H's in a graph G. Then, there is an 0{n + a{G)m + t{G)) 
time algorithm counting the number of induced H's of G for every graph H on four vertices. 

Proof. Let k, d, s, p, q, and y denote the number of induced K^s, diamonds, squares, P4's, 
paws, and claws in G, respectively Similarly, let k, d, s, q and y be the number of induced 
complements of K^^s, diamonds, squares, paws, and claws in G, respectively. Define rn = 
(2) ~ "^' ^(^) = ^ ~ d{v) — 1 for u G V{G), and 5{v, w) = d{v) — d{vw) for vw G E{G). That 
is, m is the number of edges of G, d{v) is the degree of v in G, and 5{v, w) is the number of 
vertices that are adjacent to v and not w. Finally, let S be the set obtained after executing 
the algorithm C4 in [3j. In 5, each element is a triple {v,w,L) where v,w € ViG) and L is 
a set of vertices. Then, G fulfills the following system of linear equations. 

E 0^^^)=3k + d + s (1) 

j: ("(;"■')= 6. +d (2) 

y S{v,w)6{w,v) = 4s + p (3) 

vweE{G) 

v^ fd(v) + d(w) - divw) - 2\ , , , 

Y^ (^' ^2 j = 6A; + 5d + 4s + p + 3g + 3y (5) 

J]] (i'(t;)(n-3) = 12A; + 6d + 3g + 3y (6) 

^ t)ey(G) ^ ^ 
7) " 2^ ('^2''^) =s+P + 2>k + 2d + 27s + q (9) 

?2 \ — — 

\=k + d+s+p + q + y + k + d + s + q + y (10) 

Equations ©-([I]) correspond to the first, third and fifth equations of [lO], respectively. For 
the translation between them, observe that if A is the adjacency matrix of G and G is the 
adjacency matrix of the complement of G, then A^^ = d{vw) and AG^^w = 5{v, w). Thus, G 
satisfies these equations [TO] . 

Each triple of S represents a set of non induced cycles of G in the following way. Let 
f 1, . . . , Un be an ordering of V{G) in non increasing order of degree. Then {vi, Vj, L) is a triple 
of S if and only if i < min{j, fc,/}, L C N{vi) Pi N{vj) and \L\ > 2. Therefore, (2 ) counts 
the number of non induced 4-length cycles that contain Vi and Vj and such that ViVj is not an 



edge of the cycle. Such cycles count three times each K4, and once each diamond and each 
square. Thus, ([T]) is fulfilled by G. 

The remaining equations follow analogously, by observing that: ([5]) counts, for each edge 
vw, the number of pairs of vertices that are adjacent to at least one of v and w; (0) counts, 
for each vertex v, the number of triangles of v plus one vertex; d?]) counts those graphs on 
four vertices where v has degree at least 2; ([8]) and Q count the number of pair of disjoint 
edges on G and G, respectively; and (jlOp counts the number of induced graphs on 4 vertices 
of G. 

These 10 equations are linearly independent, and by fixing the number of subgraphs H we 
can solve the whole system. As for the time complexity of computing the constant terms of 
the system, algorithm C4 takes 0{n + a{G)m) time [3J; d' = \N'\ is found in 0{n + a{G)m) 
time; the degrees of the edges and 5 are found in 0{n + a{G)m) by traversing all edge- 
neighborhoods; and all the terms in ([7|)- (|1U|) are computed in 0{n + m) time. D 

We now consider a slight modification of the 4-subgraph counting problem. Given a 
vertex v, the goal is to count the number of graphs on 4 vertices that contain v. We focus 
our attention on four types of connected graphs: K^'s, diamonds, paws, and claws. For 
i G {1,2,3}, define ki{v), di{v), qi{v), and yi{v) as the number of -ff4's, diamonds, paws, and 
claws that contain a given vertex v, where the degree of v in such an induced subgraph is i. 
In the previous section we saw that k3{v) can be computed in 0{d{v)h{G)a{G)) time. The 
following theorem shows that we can compute di{v), qi{v), and yi{v), once k3{v) is given. 

Theorem 10. There is dynamic graph data structure, consuming 0{n + m) space, with the 
following properties: 

• Both the insertion and the removal of a vertex v take 0{d{v)h{G)) time. 

• The time required for inserting the vertices of G, one at a time, into an initially empty 
instance of the data structure is 0{n + a{G)m). 

• Given k3{v), the values of di{v), qi{v), and yi{v) can he found in 0{d{v)h{G)) time, for 
every v S V{G) and i G {1, 2,3}. 

• If k^{v) is given for every v G V{G), then the number of diamonds, paws, and claws 
can be found in 0{n + a{G)m) time. 

Proof. Fix a vertex v and let ki = ki{v), di = di{v), qi = qi{v), and yi = yi{v), for i G {1, 2, 3}. 
Define 6 as in Theorem [9l Then, v fulfills the following system of linear equations. 

d3 = E ("^^7^) -3^3 d, = Y. id{wz)-l)-3ks 

2^3 = y^ d{vw)6{v,w) — 2ds 52 = /^ d{vw)6{w,v) — 2d2 

w£N{v) weN{v) 

q^ = Y^ [d'{w) - d{vw)) - 3ks - 2^2 

w£N{v) 

''5{v,w)\ \--^ f5{w,v] 



3.3 = E (**T') 



qs yi = z^ 1 2 ' ~ ^^ 



w£N{v) ^ ^ ■w€N{v) 

These equations are similar to those in Theorem |9l e.g., ( ^^'^ ) counts the number of pairs of 



vertices x,z that are both adjacent to v and vu, thus ^^weNiv) ( 2^" ) counts three times each 
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K4 that contains v, and one time each diamond that contains v with degree 3. The other 
equations fohow analogously. 

The dynamic data structure is just the /i-graph data structure with the addition that d{wz) 
is stored in data(w,z), for every edge wz. This value can be easily updated in 0{d{v)h{G)) 
time when a vertex v is either inserted or removed. Indeed, according to whether v is inserted 
or removed, we can increase or decrease in 1 the value of d{wz) for every wz S N'{v) with a 
single traversal of N'{v). On the other hand, the value of d{vw) is simply the degree of w in 
G[N(v)], for every w £ N{v). Then, the values of di{v), qi{v), and yi{v) can be obtained in 
0{d{v)h{G)) time by solving the system of equations above, once k^{v) is given. Furthermore, 
by solving the above equations for all the vertices, we can compute the number of diamonds, 
paws, and claws can be found in 0{n + a{G)m) time, by Corollary [5j D 

5.3 Dynamic Recognition of Diamond-Free Graphs 

Recall that the graph that is obtained by removing one edge from a complete graph of four 
vertices is called a diamond. A diamond-free graph is a graph that contains no induced 
diamond. Diamond-free graphs appear in many contexts; for example in the study of perfect 
graphs HEIIIB]- 

In [To], Kloks et al. showed how to find an induced diamond in Oivr?''^ + n'^) time, if one 
exists. The fast matrix multiplication algorithm is used in one of the steps of this algorithm, 
which explains why n^ is a term of the complexity order. However, the fast matrix multipli- 
cation can be avoided while improving the time complexity to Oirr?''^^ time, as it is shown 
by Eisenbrand and Grandoni [6]. Talmaciu and Nechita [T^ devised a recognition algorithm 
based on decompositions, but they claim that in the worst case the time required by their 
algorithm is not better than the one by Kloks et al. Note that Theorem [9] implies that there 
is an 0(a(G)^m) time algorithm for recognizing whether a graph is a diamond-free graph, 
improving over the previous algorithms for some sparse graphs. Finally, Vassilevska [l_9j used 
the algorithm by Eisenbrand and Grandoni to find an induced K]. \ e in a graph. A K^^ \ e is a 
complete graph on k vertices, minus one edge. For every even A;, the algorithm by Vassilevska 
takes O [d{n, m)m^ ''"^^ time, where d{n, m) is the time required to find a diamond in a 
graph with n vertices and m edges. Thus, this algorithm is implicitly improved with each 
improvement on d{n,m). 

The algorithm by Kloks et al. is based on the fact that a graph G is diamond-free if and 
only if G[N{v)] is a disjoint union of maximal cliques, for every v G V{G). Testing whether 
a graph is a disjoint union of cliques takes linear time, and we saw in Section [5] how to 
compute the family {G[A^(w)]}^gv'{G) ™ 0{a{G)m) time. Therefore, by using the /i-graph 
data structure, the algorithm by Kloks et al. can be implemented so as to run in 0{a{G)m) 
time, improving over the algorithm by Eisenbrand and Grandoni and the algorithm implied 
by Theorem [9j 

In this section we develop a dynamic data structure for maintaining diamond-free graphs. 
This data structure can also be used to find an induced diamond of a static graph G in 
0{a{G)7n) time, if one exists. As a by-product, the data structure can be used to query the 
maximal cliques of the dynamic diamond-free graph in constant time. The data structure is 
based on this well known theorem about diamond-free graphs. 

Theorem 11. A graph is a diamond-free graph if and only if every edge belongs to exactly 
one maximal clique. 



11 



An /i-graph data structure is used to represent a diamond-free graph G. Also, the family 
of non singleton cliques of G is stored in the dynamic diamond-free data structure. Each non 
singleton clique is in turn implemented as the set of edges that belong to the clique. Also, the 
pointer data(t;tt;) references the unique clique G^w of the family that contains vw, for every 
vw G E{G). Finally, each clique G is associated with a counter c{G) that is initialized to 0. 
The purpose of c{G) is to count the number of neighbors of v inside C, when a vertex v is 
inserted into G. 

We are now ready to discuss the operations allowed by the data structure. Suppose first 
that G is a diamond-free graph and v V{G) is to be inserted into G. Vertex v is given 
with its neighborhood set N{v), which also defines its edge-neighborhood N'{v). Say that v 
is edge-adjacent to a clique C of G when E{G) n N'{v) / 0, while v is fully edge-adjacent to 
G when E{G) n N'{v) = E{G). The following theorem shows how to insert v into G. 

Theorem 12. The graph GU {v} is diamond-free if and only if the following two statements 
hold for every maximal clique G to which v is edge- adjacent. 

1. V is fully edge-adjacent to G , and 

2. if V is edge-adjacent to a maximal clique G' ^ G , then V{G') n V{G) = 0. 

Proof. If V is not fully edge- adjacent to C, then there is some vertex u £ G which is not 
adjacent to v. Since v is edge-adjacent to C, then there is an edge wz G N'{v) n E{G). But 
then, u, V, w, z induce a diamond in G. Suppose now that v is edge-adjacent to a maximal 
clique G' ^ G that contains some vertex u G G. Since G and C" are maximal cliques, it 
follows that there is some vertex w G G which is not adjacent to a vertex z £ G' . Thus, 
u,v,w,z induce a diamond. 

For the converse, suppose that G U {v} is not a diamond-free graph. Since G is diamond- 
free, then there are three vertices u, w, z such that together with v induce a diamond in G. 
If u and V are not adjacent, then u,w,z belong to some maximal clique G of G. Since v is 
adjacent to w, z but not to u, we obtain that v is edge-adjacent but not fully edge-adjacent 
to G. So, we may assume that v is adjacent to u, w, z, and that w and z are not adjacent. 
But then, uw and uz belong to different maximal cliques G and C". Thus, v is edge-adjacent 
to G and G' and u G V{G) n V{G'). D 

Algorithm [21 which is obtained from Theorem 1 12[ can be used to decide whether G U {v} 
is a diamond-free graph. 

Algorithm 2. Insertion of a vertex v into a diamond-free graph G 

1. Remove the mark from every w G V{G). 

2. For each maximal clique C of G to which v is edge-adjacent: 

3. If V is not fully edge-adjacent to C, then output "G U {v} is not diamond-free." 

4. For each maximal clique G of G to which v is edge-adjacent: 

5. If there is some marked w £ G, then output "G U {v} is not diamond-free." 
Otherwise, mark every w £ G. 

6. Output "G U {v} is a diamond-free graph" . 
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Observe that after Gu{v} is claimed to be a diamond-free graph in Step[6l all the vertices 
of every edge- adjacent maximal clique C have a mark, and every vertex w & C was traversed 
and marked only once. Then, v is fully-adjacent to C and no vertex of C belongs to other 
clique to which v is edge-adjacent. Therefore, by Theorem 1121 the algorithm is correct. 

Discuss the implementation of Algorithm [2j The input of the algorithm is formed by the 
graph G, the vertex v, and the set N{v) of neighbors of v in G. The first step is to insert v 
into the /i-graph data structure of G in 0{d{v)h{G)) time. We consider that w is marked if 
and only if data(i;t(;) 7^ NULL. Hence, Step [1] of the algorithm is executed in 0{d{v)) time, 
by traversing once the set N{v) to initialize data(fty). 

Before traversing each edge-adjacent maximal clique in Steps [2] and HI we compute N'{v), 
in 0{d{v)h{G)) time, as in Section[H The family C = {Cwz}wz(^N'(v) is computed by iteratively 
inserting Cyjz into C, for every wz G N'{v). While C is generated, \N'{v) n E{Cyjz)\ can be 
computed by increasing c{Gyjz) by 1 when wz is first traversed. This operation takes 0(1) 
time per wz E N'(v), thus it takes 0{d{v)h[G)) total time. Once C is computed, each maximal 
clique C G C is traversed. 

If c(C) = \N' {v)f\E{C)\ / \E{C)\, then v is not fully edge-adjacent to C and the algorithm 
stops with a failure message in Step [3l Otherwise, every maximal clique C ^ C has to be 
traversed once again to update da.tz.{vw) so as to reference C„^, for every w £ G. Instead of 
doing this, we traverse N'{v) and, for every edge wz, we set data(z;i(;) and data(t'2;) to point 
to Cwz, because v is fully adjacent to C^z- If when updating C^w (resp. C^z) to point to C we 
discover that Gyw (resp. C^z) points to C" with G / C", then the algorithm stops in failure 
as in StepO If C^z is maintained together with wz in N'[v), then each edge gets traversed 
in 0(1) time. Thus, the time complexity of both loops is 0{d'{v)) = 0{d{v)h{G)). 

After the main loops are completed, the algorithm claims that G U {u} is a diamond-free 
graph. Before another vertex can be inserted, the data structure that represents G has to be 
updated into a data structure that represents GU {v}. The set of edges {vw \ w G N{v)} 
can be split into two types. Those whose pointer data references some maximal clique of G 
and those whose data is still NULL. For those edges that reference some maximal clique, 
the edge has to be inserted into Gyyj. For those edges vw with data = NULL, a new clique 
has to be created that contains only the edge vw. This clique has to be inserted to the set 
of maximal cliques of G U {u } and data(t'u;) has to be updated accordingly. All this can 
be done in 0{d{v)) time. Finally, c{G) has to be updated to for every C G C, again in 
0{d'{v)) = 0{d{v)h{G)) time. Therefore, the insertion of a new vertex takes 0{d{v)h{G)) 
time as desired. 

Observe that if G U {v} is a diamond-free graph, then its family of maximal cliques is 
obtained in 0(1) time and the maximal clique to which an edge vw belongs can be queried 
in 0(1) time. Algorithm [2] can also be modified so as to output an induced diamond in 
0{d{v)h{G)) time, when G U {v} is not diamond-free. Consider the two alternatives for the 
algorithm to stop in failure. First, if v is not fully edge-adjacent to C G C, then there is some 
vertex w £ G that is not adjacent to v. In this case, v,w and the endpoints of an edge in 
N'{v) n G induce a diamond. To find w, traverse the edges of G and query if each endpoint is 
adjacent to v. (For this step, we may use the marks used in the computation of N'{v), or we 
may re-compute them with a single traversal of N{v).) The first vertex that is not adjacent 
to V is taken as w. To find an edge in N'{v) n G, traverse N'{v) until some edge of G is 
reached. Thus, the certificate in this case can be found in 0(d'{v)) time. The second reason 
for the algorithm to stop in failure is that v is fully edge-adjacent to G and C", and vw is 
marked with G' while trying to mark it with G ^ G'. In this case, we ought to find two non 
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adjacent vertices of CU C". As in the proof of Theorem 1121 these two vertices together with 
V and w induce a diamond. To find these two vertices observe that w is the unique vertex in 
C n C" and that all the vertices of C \ {w} are not adjacent to the vertices of C \ {w}. Thus, 
the two non adjacent vertices of C U C" can be found in 0(1) time by traversing at most one 
edge of both C and C". 

The dynamic data structure supports also the removal of vertices and the insertion and 
removal of edges. For the removal of a vertex v, note that G \ {v} is always a diamond-free 
graph. Thus, all we need to do is to remove v from the /i-graph data structure, and to remove 
vw from Cjjyi], for every w € N{v). As explained in Sectional the former operation takes 
0{d{v)h{G)) time, while the latter takes 0{d{v)) time. With respect to the insertion of an 
edge vw, graph G U {vw} is diamond-free if and only if \N{v) n N{w)\ < 1 and, if there is 
some z £ N{v) Ci N{w), then dcivz) = dciwz) = 0. We can compute \N{v) n N{w)\ in 
0{d{v) + d{w)) time, while, for z G N{v) n N{w), we can access data(uz) and data(it;z) 
in 0(1) time so as to see whether these maximal cliques have exactly one edge. If so, the 
maximal cliques C^z and C^z need to be merged, again in 0(1) time. Thus, the insertion of 
vw takes 0{d{v) + d{w)) time. Finally, for vw £ E{G), the graph G\{vw} is diamond-free if 
and only if d{vw) < 1, i.e., if C^w has at most three edges. If d{vw) = 1, then C^w has to be 
split into two maximal cliques of O \ {vw} in 0(1) time. Therefore, the removal of an edge 
takes 0{m.m{di{v),di{w)} + hi{v) + hi{w)) because we need to remove vw from the /i-graph 
data structure. 

Finally, the diamond-free data structure can be used to test whether a graph G is diamond- 
free, just by iteratively inserting the vertices of G into the data structure. At each step, the 
operations of greater complexity are the insertion of the new vertex into the /i-graph data 
structure, and the computation of N'{v). As we have discussed in Section UJ the total time 
cost for these operations is 0{a{G)m). Thus, this algorithm runs as fast as the improvement 
of the algorithm by Kloks et al. discussed before. 

5.4 Simple, Simplicial, and Dominated Vertices of Dynamic Graphs 

A vertex v is dominated by a vertex w if N[v] C N[w]. Equivalently, v is dominated by 
w G N{v) if d{v) — d{vw) = 1. We say that v and w are comparable if either v is dominated 
hy w or w is dominated by v. If v is dominated by all its neighbors, then u is a smplicial 
vertex, while if t> is a simplicial vertex and every pair of neighbors are comparable, then v is 
a simple vertex. 

In [10], Kloks et al. showed how to compute the set of simplicial vertices in 0{n + 
j^2aj/(aj+i)'j _ (9 (t^ _j_ 777,1-41^ time, usiug the fast matrix multiplication algorithm. With the 
/i-graph data structure, we can find all the simple, simplicial, and dominated vertices in 
0(a(G)m,) time, as follows. First, find the degree d{vw) for every vw G E{G) in 0{n+a{G)m) 
time, as discussed in Section HI Next, for each vertex f , find the set of vertices D[v) that 
dominate v, by testing whether d{v) — d{vw) = 1, for every w G N{v). Clearly, if D{v) ^ 0, 
then t> is a dominated vertex, while if |-D(t>)| = d{v), then t> is a simplicial vertex. To de- 
termine if a simplicial vertex v is simple, it is enough to check whether z G D{w), for each 
edge wz G N'{v) with d{w) < d{z). As discussed in Section [H we can traverse all the edge- 
neighborhoods N'{v) in 0{a{G)m) time, thus simple, simplicial, and dominated vertices can 
be found in 0{a{G)m) time. 

In this section, we show how can the /i-graph data structure be used to maintain the 
simple, simplicial, and dominated vertices, while vertices are inserted to or removed from a 
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dynamic graph. Let G be a graph implemented with the /i-graph data structure, where d{wz) 
is stored in data(i(;, z) for every edge wz, and let D be the family of dominated vertices. 
Suppose that a new vertex v with neighborhood N[v) C V{G) is to be inserted into G, and 
that we want to update D so as to store the dominated vertices of GU {u}. The next lemma 
shows how to find the new dominated vertices. 

Lemma 13. A vertex w ^ v is dominated in H = GU {v} if and only if one of the following 
statements is true: 

• w ^ N{v) and w is dominated in G, 

• dniw) — dnivw) = 1, or 

• dniw) — dniwz) = 1, for some wz G N'fj{v). 

The first step is to insert v into G and to compute d{vw) for every w G N{v). Both steps 
take 0{d{v)h{G)) time, as discussed in Section[l]and Theorem [TOl Next, we update the set D. 
By the lemma above, we need not consider the vertices outside N{v). To find those w G N{v) 
that are dominated by v, we traverse N{v) while checking whether d{w) — d{vw) = 1. Next, we 
find all the other dominated vertices, by checking the values of d{w) — d{wz) and d{z) — d{wz), 
for every wz G N'{v). Finally, we remove from D those neighbors of v that are not longer 
dominated, and we insert v if there is some w such that d{v) — d{vw) = 1. Since N'{v) is 
computed in 0{d{v)h{G)) time, the whole procedure takes 0{d{v)h{G)) time. 

A similar procedure can be used to update the family S of simplicial vertices of G when 
V is inserted into G. In this case, the simplicial vertices are found as in the next lemma. 

Lemma 14. A vertex w ^ v is simplicial in H = GU {v} if and only if one of the following 
statements is true: 

• w ^ N{v) and w is simplicial in G, or 

• w is simplicial in G, and d{w) — d{vw) = 1. 

Again, begin with the insertion of v into G and the computation of d{vw) for every 
w G N{v), in 0[d{v)h{G)) time. In this case, we update S by first traversing each w G N{y) 
and checking whether w £ S and d{w) — d{vw) = 1. On the other hand, we insert v into S if 
and only if d{v) — d{vw) = 1 for every w G N{v). The time required by these operations is 
0{d{v)), once v was inserted into G. Thus, the update of S takes 0{d{v)h{G)) time. 

Updating the family Q of simple vertices is not as simple as updating the families D and 
S. The reason is that we can no longer skip the vertices outside N{v). To provide an efficient 
update of Q, we store in data(a;), for each vertex x, the number fi{x) of edges wz G N'(x) such 
that w and z are not comparable. So, x is simple if and only if x is simplicial and /x(a;) = 0. 
We can find out the value of fj,{v) in 0{d{v)h{G)) time, by traversing every wz G N'{v) and 
checking whether ra.m.{d{w) , d{z)} — d{wz) = 1. The update of all the other values of /x is 
based on the following lemma. 

Lemma 15. Two vertices w and z of G are comparable in G U {v} if and only if they are 
comparable in G and either {w, z} C N{v) or {w, z} D N{v) = 0. 

The update of /i(x) for every x 7^ v is done as in Algorithm [3j 
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Algorithm 3. Update of /x after the msertion of v. 

1. For each w € N{v): 

2. For each z G H{w): 

3. If z ^ N{v), dciw) — dciwz) = 1 and dciz) — dciwz) > 1, then: 

4. {w and z were comparable before inserting v, but now they are not} 

5. Set fj,{x) = fi{x) + 1, for every x G N{wz). 

Algorithm [3] requires 0{m) time in the worst case. However, suppose that Algorithm [3] is 
iteratively executed for computing the simple vertices as in Algorithm|H In such an algorithm, 
the condition of the inner loop of Algorithm [3] is executed at most once for each edge wz. 
Indeed, if w and z are not comparable prior the insertion of u , then they are not comparable 
after the insertion of v. On the other hand, when wz meets the condition of the inner loop, we 
know that dciw) < dciz). Hence, as discussed in Section HI the time required by Algorithm [J] 
is 0{n + a{G)m) for the update of the /i-graph data structure and d{vw), and 

^ ""+ X] 'am\{d{v),d{z)]\ =0{n + a{G)m) 

\ wzeE{G) J 

for the update of /i. 

Algorithm 4. Iterative update of ^ for a graph G. 

1. Let vi, . . . ,Vn he an ordering of V{G), and G' be an empty graph. 

2. For i = 1, . . . ,n: insert Vi into G' while executing Algorithm [3l 

As a corollary, we obtain that the simple, simplicial, and dominated vertices of G can be 
found in 0{n + a{G)m) time. The process for updating the sets of simple, simplicial, and 
dominated vertices when a vertex is removed is similar. When v is removed, again we should 
not consider those vertices outside N{v) for the update of D and S. A vertex in N{v) is 
dominated in G\ {t;} if, in G\{v}, d{w) — d{wz) = 1 for some z £ H{w), while it is simplicial 
in G \ {v} if L{w) = and d{w) — d{wz) = 1 for every z G H{w). Thus, we can update D 
and S in 0{d{v)h{G)) time. Finally, we can update jjl in 0{m) time when v is removed by 
applying Algorithm [Sj which is the inverse of Algorithm [3l 

Algorithm 5. Update of ^ after the insertion of v. 

1. For each w G N{v): 

2. For each z G H{w): 

3. If z ^ A(t>), dciw) — dciwz) = 2 and dciz) — dciwz) > 1, then: 

4. {w and z were not comparable before inserting v, but now they are} 

5. Set fi{x) = fi{x) — 1, for every x G N{wz). 
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The removal operations can be used to improve the best known algorithm for the recog- 
nition of cop-win and strongly perfect graphs, for sparse graphs. 

A cop-win order of a graph G is an ordering vi, . . . ,Vn of V{G) such that Vi is dominated 
in the subgraph induced by ijj, . . . , z;„, for 1 <i <n. A graph that admits a cop- win order is 
a cop-win graph. The cop- win name comes from the fact that cop-win graphs are precisely the 
graphs in which a cop can always catch the robber in a pursuit game [Jl]. This class has been 
introduced in [12], cf. [2j- Cop-win graphs are also known in the literature under the name 
of dismantlable graphs ^3j, and they are a main tool in the study of clique graphs [16] ■ The 
currently best algorithms for recognizing cop-win graphs run in 0{nm) time or in 0{n^ / log n) 
time [15]. 

A dismantling of a graph G is a graph H obtained by iteratively removing one dominated 
vertex of G, until no more dominated vertices remain. It is not hard to see that all the 
dismantlings of G are isomorphic. Using the /i-graph data structure, we can compute the 
dismantling of a graph, and the cop-win order of a cop-win graph, in 0{a{G)m) time easily. 
First, find the set D of dominated vertices in 0[n + a{G)m) time (cf. above). Then, while 
D 7^ 0, choose a vertex D, and remove it from G while updating D as explained above. The 
graph obtained after this procedure is the dismantling H of G. If H has a unique vertex, 
then G has a cop-win order, given by the order in which the vertices were removed by the 
algorithm. This algorithm takes 0{n -\- a{G)m) time, as discussed above. 

A simple elimination ordering of a graph G is an ordering vi, . . . ,Vn of V{G) such that Vi 
is a simple vertex in the subgraph of G induced by Uj, . . . , v„. The family of graphs that admit 
a simple elimination ordering is precisely the family of strongly chordal graphs [^j . Strongly 
chordal graphs were introduced as a subclass of chordal graph for which the domination 
problem is solvable in polynomial time [8J. The best known algorithms for computing a simple 
elimination ordering run in either 0{n'^) or 0{{n-\-m) logn) time |14j . These algorithms work 
by finding a doubly lexical ordering of the adjacency matrix of the graph, and then testing 
if this sorted matrix contains some forbidden structure. Our approach, instead, is based on 
iteratively finding a simple elimination ordering by iteratively removing the simple vertices. 

Every strongly chordal graph G has at least one simple vertex, and G \ {v} is strongly 
chordal for every v G y{G) [8j. Thus, we can compute a simple elimination ordering of G 
by removing the simple vertices in any order. That is, first find the set Q of simple vertices 
in 0{n -\- a{G)m) time (cf. above). Then, while Q 7^ 0, choose a vertex in Q, and remove it 
from G while updating Q with Algorithm [5j If all the vertices are removed from Q, then the 
removal order given by the algorithm is a simple elimination ordering of G. As discussed for 
Algorithm HI the inner loop of Algorithm \5\ is executed at most once for each edge (cf . above) . 
Therefore, the simple elimination ordering is computed in 0{n-\-a{G)m) time, improving the 
previous best algorithms for graphs with low arboricity. 
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